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Recently, attempts have been made to understand the ap- 
parent near coincidence of the present dark energy and matter 
energy in terms of a dynamical attractor-like solution for the 
evolution of a "quintessence" scalar field. In these models the 
field couples with the dominant constituent and only acts like 
a cosmological constant after the onset of the matter domi- 
nated epoch. A generic feature of such solutions, however, is 
the possibility of significant energy density in the scalar field 
during the radiation dominated epoch. This possibility is even 
greater if the quintessence field begins in a kinetic-dominated 
regime generated at the end of "quintessential infiation." As 
such, these models can affect, and therefore be constrained by, 
primordial nucleosynthesis and the epoch of photon decou- 
pling. Here we analyze one popular form for the quintessence 
field (with and without a supergravity correction) and quan- 
tify constraints on the allowed initial conditions and param- 
eters for the effective potential. We also deduce constraints 
on the epoch of matter creation at the end of quintessential 
infiation. 

PACS numbers; 98.80.Cq, 98.65.Dx, 98.70.Vc 



I. INTRODUCTION 

Recent observations [^^J of Type la supernovae at in- 
termediate redshift, along with complementary observa- 
tional constraints at low and intermediate redshift, as 
well as the power spectrum of the cosmic microwave back- 
ground all indicate |^ that the universe may be presently 
accelerating due to the influence of a dominant dark en- 
ergy with a negative pressure. The simplest interpreta- 
tion of this dark energy is a cosmological constant for 
which the equation of state is u! = P/ p — —l. A second 
possibility is derived from the so-called "quintessence" 
models. In this case the dark energy is the result of a 
scalar field Q slowly evolving along an effective potential 
V{Q). The equation of state is negative — 1 < wq < 0, 
but not necessarily constant. 

Introducing either of these paradigms, however, leads 
inevitably to two nagging questions. One is a fine tuning 
problem as to why the present dark energy is so small 
compared to the natural scales of high-energy physics. 



The second is a cosmic coincidence problem as to why 
the universe has conspired to produce nearly equivalent 
energy content in matter and dark energy at the present 
time. 

Attempts have been made P-|l0|] to reformulate 
this quandary by introducing specific forms of the 
quintessence effective potential whereby a tracker field 
Q evolves according to an attractor-like solution to the 
equations of motion. That is, for a wide variety of ini- 
tial conditions, the solutions for Q and Q rapidly ap- 
proach a common evolutionary track. The nice feature 
of these solutions is that they lead naturally to a cross 
over from an earlier radiation-dominated solution to one 
in which the universe is dominated by a small dark energy 
at late times. Another interesting possible feature is that 
such models might also naturally arise during matter cre- 
ation at the end of an earlier "quintessential" inflationary 
epoch . In this case, the Q field emerges in a kinetic- 
dominated regime at energy densities above the tracker 
solution. 

It is not yet clear, however, that these mod- 

els have altogether solved the fine-tuning and cosmic- 
coincidence problems. Nevertheless, several such tracker 
fields have been proposed jl^. Although there is some 
difficulty in aligning quintessence models and string the- 
ory , the form for the effective potentials may at least 
be suggested by particle physics models with dynamical 
symmetry breaking, by nonperturbative effects by 
generic kinetic terms "fc-essence" in an effective action 
describing the moduli and massless degrees of freedom in 
string and supergravity theories H^H], or by static and 
moving branes in a dilaton gravity background 

A general feature of all such solutions, however, is that 
at least the possibility exists for a significant contribu- 
tion of the Q field to the total energy density during 
the radiation-dominated or photon decoupling epochs as 
well as the present matter-dominated epoch. The yields 
of primordial nucleosynthesis and the power spectrum of 
the CMB are strongly affected by the background en- 
ergy density and universal expansion rate. Therefore, 
the possibility exists to utilize primordial nucleosynthe- 
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sis and the CMB power spectrum to constrain otherwise 
viable quintessence or fc-essence models. 

Observational constraints on such quintessence mod- 
els have been of considerable recent interest. For exam- 
ple, in jl^ the constraints on an M-theory motivated 
exponential form with a quadratic prefactor for the 
quintessence effective potential were considered. In this 
model, the Q field closely follows the background field. 
In a study was made of the effect on primordial nu- 
cleosynthesis of extended quintessence with nonminimal 
coupling. In such models the quintessence field couples 
with the scalar curvature. It can then act to slow the 
expansion and affect the yields of primordial nucleosyn- 
thesis. In the present work, we consider the power-law 
effective potential and its supergravity-corrected form as 
described below. We construct a detailed mapping of 
the allowed parameter space for these quintessence mod- 
els. We also consider constraints these considerations 
place on the epoch of matter generation at the end of 
quintessential inflation. 



Kahler potential yields an extra factor of exp{'iQ'^ / 2rn^i} 
pof . This comes about by imposing the condition that 
the expectation value of the superpotential vanishes. The 
Ratra potential thus becomes 
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(Q) ^ m(4+")q-" X exp{3Qy2mli) , (4) 



where the exponential correction becomes largest near 
the present time as Q — > rripi. This supergravity mo- 
tivated effective potential is known as the SUGRA po- 
tential. The fact that this potential has a minimum 
for Q = ^JotlZm^i changes the dynamics. It causes the 
present value of wq to evolve to a cosmological constant 
(wq ~ — 1) much quicker than for the bare power-law 
potential 

The quintessence field Q obeys the equation of motion 



Q + 3HQ + dV/dQ ^ 



(5) 



where the Hubble parameter H is given from the solution 
to the Friedmann equation, 



II. QUINTESSENCE FIELD 

A variety of quintessence effective potentials or 
/c-essence effective actions (g-gll can be found in the lit- 
erature. In this paper, however, we will not consider 
/c-essence effective actions as we have found that the 
attractor-like solution is limited to a prohibitively small 
range of parameters of the effective Lagrangian. We will, 
however, consider kinetic dominated quintessence mod- 
els. Here, we concentrate on what is a frequently invoked 
form for the effective potential of the tracker field, i.e. an 
inverse power law such as originally analyzed by Ratra 
and Peebles ITgl, 



(1) 



where, M and a are parameters. The parameter M in 
these potentials is fixed by the condition that 17q = 0.7 
at present. Therefore, 



Pq(0) = 0.7/9c(0) = 5.7h^ X 10"^' GeV^ (2) 



and 



M« Pq(0)Q" 



l/{a+i) 



(3) 



If Q is presently near the Planck mass and a is too small 
(say a^2), this implies a small value Q for M which in 
a sense reintroduces the fine tuning problem. 

We will also consider a modified form of V{Q) as 
proposed by ||2^ based upon the condition that the 
quintessence fields be part of supergravity models. The 
rewriting of the effective potential in supergravity de- 
pends upon choices of the Kahler potential . The flat 



-\-{pb+Pq) , 



(6) 



where rUpi = {8ttG/3)-^^^ = 4.2 x 10^*^ GeV, pg is the 
energy density in background radiation and matter, and 
a is the cosmic scale factor. 

For the simple inverse power-law potential, it can be 
shown that the tracker solution maintains the con- 
dition 



d^V/dQ^ = (9/2)(l - + l)/a]H^ 



(7) 



As the Q field evolves, its contribution to the energy 
density is given by 



PQ = ^ + V{Q) , 



(8) 



which may or may not be comparable to the energy den- 
sity in radiation during the nucleosynthesis or photon- 
decoupling epoch depending upon the parameters and 
initial conditions. 

The quintessence initial conditions are probably set in 
place near the inflation epoch. By the time of the big 
bang nucleosynthesis (BBN) epoch, many of the possible 
initial conditions will have already achieved the tracker 
solution. However, for initial conditions sufficiently re- 
moved from the tracker solution, it is quite possible that 
the tracker solution has not yet been achieved by the time 
of BBN. Such possibilities are illustrated schematically in 
Figure 0. These correspond to cases in which the initial 
energy density falls above (curve A), near (curve B) or 
below (curve C) the tracker solution. 
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III. NUCLEOSYNTHESIS CONSTRAINT 
A. Quintessence and BBN 

There are several paradigms of possible interest for 
constraint by primordial nucleosynthesis. These depend 
upon the initial values for the energy density in the Q 
field. In any of these cases, the energy density in pq can 
be constrained by the ratio pq / ps during the BBN epoch 
at 0.01<T<1 MeV, 10^ <z as shown in Figure |l]. 

If the initial conditions are sufficiently close to the 
tracker solution pq before nucleosynthesis, then pq ^ pq 
during nucleosynthesis and the tracker solution speci- 
fies the energy density in the Q field. This situation 
is similar to curve B in Figure 1. This may be the 
most likely scenario. Along the tracker solution, pq di- 
minishes in a slightly different way than the radiation- 
dominant background energy density. For example, as 
long as Pq « pb, the Q- field decays as 

PQ (X a-3(i+'"«) , (9) 

with 

wq = {awB-2)/{a + 2)<WB ■ (10) 

The equation of state wq is only equal to the back- 
ground equation of state wb in the limit a ^ oo. This is 
equivalent to the exponential potential. Nevertheless, the 
tracker solution does not deviate much from ps, even at 
high redshift for most values of a considered here. Hence, 
one can characterize the nucleosynthesis results by the 
(nearly constant) ratio pq/pB during the BBN epoch. 

If the energy density in the tracker solution is close 
to the background energy density, the nucleosynthesis 
will be affected by the increased expansion rate from the 
increased total energy density. Such a situation occurs 
for large values of the power-law exponent a. 

A second possibility is that the energy density pq 
could exceed the tracker solution and be comparable to 
or greater than the background energy density during 
primordial nucleosynthesis. This situation is something 
like curve A in Figure 0. The kinetic energy in the Q 
field dominates over the potential energy contribution to 
Pq and wq = -1-1 so that the kinetic energy density di- 
minishes as a~^. In this case there could be a significant 
contribution from pq during nucleosynthesis as the Q 
field approaches the tracker solution. The strongest con- 
straints on this case would arise when pq is comparable 
to the background energy density near the time of the 
weak-reaction freese out, while the later nuclear- reaction 
epoch might be unaffected. This case A is particularly 
interesting as this kinetic-dominated evolution could be 
generated by an earlier quintessential inflation epoch | [Tl| ] 
as described below. 

A final possibility might occur if the Q field approaches 
the tracker solution from below as indicated by curve C 



in Figure [l|. In this case, the tracker solution may be 
achieved after the the BBN epoch so that a small pq 
during BBN is easily guaranteed. In such models, how- 
ever, the ultimate tracker curve might have a large energy 
density at a later epoch which could be constrained by 
the later CMB epoch as described below. 

B. Light-Element Constraints 

The primordial light-element abundances deduced 
from observations have been reviewed by a number of 
recent papers |pl|-^4| . There are several outstanding un- 
certainties. For primordial helium there is an uncertainty 
due to the fact that the deduced abundances tend to re- 
side in one of two possible values. Hence, for ^He we 
adopt a rather conservative range: 

0.226 <Yp< 0.247. 

For deuterium there is a similar possibility for either a 
high or low value. Here, however, we adopt the generally 
accepted [^|3| low values of 

2.9 X 10"^ < D/H < 4.0 X 10"^ 

The primordial lithium abundance can be inferred from 
the surface abundances of old halo stars. There is, 
however, some uncertainty from the possibility that old 
halo stars may have gradually depleted their primordial 
lithium. Because of this we do not consider primordial 
lithium in the discussions here. The adoption of the low 
deuterium abundance forces the BBN primordial helium 
to be near its upper limit. 

Adding energy density from the Q field tends to in- 
crease the universal expansion rate. Consequently, the 
weak reaction rates freeze out at a higher temperature 
Tu,. This fixes the neutron to proton ratio (n/p ~ 
exp[(mp — m„)/r^]) at a larger value. Since most of the 
free neutrons are converted into *He, the primordial he- 
lium production is increased. Also, since the epoch of 
nuclear reactions is shortened, the efficiency of burning 
deuterium into helium is diminished and the deuterium 
abundance also increases. Hence, very little excess en- 
ergy density from the Q field is allowed. 

Figure |^ summarizes allowed values of pg/ps at T =1 
MeV (z w 10^") based upon the nucleosynthesis con- 
straints. The region to the right of the curve labeled D/H 
corresponds to models in which the primordial deuterium 
constraint is satisfied, D/H < 4.0 x 10^^. The region be- 
low the line labeled Yp corresponds to Yp < 0.247. The 
hatched region summarizes allowed values for the energy 
density in the quintessence field during the nucleosynthe- 
sis epoch. 

This figure is similar to that obtained in Freese et al. 
p5[ based upon the light-element constraints available at 
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that time. They similarly considered dark energy densi- 
ties which scale proportionally to the background radia- 
tion energy density. However, in their study a somewhat 
larger range of possible dark energies was allowed due to 
the larger uncertainties in the observed primordial abun- 
dances and the neutron half life at that time. 

In the present work we deduce an absolute upper limit 
of 5.6% of the background radiation energy density al- 
lowed in the quintessence field. This maximum contribu- 
tion is only allowed for 7710 ~ 4.75 or rif,/i^ « 0.017. 
A smaller contribution is allowed for other values of 
ryio- Indeed, this optimum 7710 value is 4cr less than 
the value implied by the cosmic deuterium abundance 
||,|| ni,h^ = 0.020 ± 0.001 (la) {r]io = 5.46 ± 0.27). 
The most recent independent determinations of fli,h^ 
from high-resolution measurements of the power spec- 
trum of fluctuations in the cosmic microwave background 
favor a value even higher. Both the BOOMERANG ^] 
and DASI Q data sets imply nth^ = 0.022ig °g^ (Ict) 
(7710 — 6.001q g5). The deuterium and CMB constraints 
together demand that 7710 > 5.19 which would limit the 
allowed contribution from the Q field to < 2% of the 
background energy density. 

The most restrictive CMB constraint on ilt,h^ de- 
rives from demanding a flat universe {Qtot — 1-0), and 
marginalizing the likelihood function over all parameters 
with assumed Gaussian priors |p6| based upon measure- 
ments of large-scale structure and Type la supernovae. 
This gives Qbh"^ = 0.023±0.003 (Ict). If one adopts this as 
a most extreme case, then r2f,/i^ > 0.020. This would cor- 
respond to ?7io > 5.46. From Figure 2 this would imply a 
much more stringent constraint that only about 0.1% of 
the background energy density could be contributed by 
the Q field. Of course, this is only a Ict constraint, and 
it is questionable as to whether the upper limit to Yp is 
well enough established to rule out a contribution to the 
energy density at the 0.1% level. For the remainder of 
this paper we will adopt the more conservative constraint 
of 5.6%. Nevertheless, it is of interest to explore how the 
quintessence parameters allowed by BBN might improve 
should the constraints from BBN ever be so tightly de- 
fined. Therefore, we will consider 0.1% as a conceivable 
limit that demonstrates the sensitivity to BBN. Even the 
most conservative 5.6% limit adopted here corresponds to 
only about half of the energy density allowed in |^ for 
3 neutrino flavors. 

IV. EQUATION OF STATE CONSTRAINT 

Another constraint on allowed parameters for the 
quintessence field derives from the simple requirement 
that the Q field behaves like dark energy during the 
present matter-dominated epoch. Stated another way, 
the equation of state should be sufficiently negative, i.e. 
Wq = Pq/pq < by the present time. 



In general, wq is a time-dependent quantity. The en- 
ergy density and pressure in the Q field can be written 
PQ = W + ■^(Q) and Pq = - V{Q). This gives, 

WQ^l~2V{Q)/pQ , (11) 

where the time dependence derives from the evolution of 
V{Q) and p{Q)- A comprehensive study of the observa- 
tional constraints on the present value of wq has been re- 
cently summarized by Wang et al. ^ . We adopt the same 
constraints deduced in that paper. Mindful of systematic 
errors, they adopted a most conservative approach based 
upon progressively less reliable data sets. Using the most 
reliable current low-redshift and microwave background 
measurements, they deduce limits of — 1 < wq < —0.2 
at the 2ct level. Factoring in the constraint from Type 
la supernovae reduces the range for the equation of state 
to — 1 < wq < —0.4. This range derives from a concor- 
dance analysis of models consistent with each observa- 
tional constraint at the 2ct level or better. A combined 
maximum likelihood analysis suggests a smaller range of 
—0.8 < wq < —0.6 for quintessence models which fol- 
low the tracker solution, though wq w — 1 is still allowed 
in models with nearly a constant dark energy. In what 
follows, we also invoke these same three possible limit 
ranges for the present value of wq . 

As we shall see, these limits place a strong constraint 
on the bare Ratra power-law potential for almost any 
value of a. However, the SUGRA corrected form of V{Q) 
is only slightly constrained. 

V. CMB CONSTRAINT 

A third constraint on the quintessence field arises from 
its effect on the epoch of photon decoupling. There are 
two effects to be considered. 

A. Look-back effect 

One is the effect of the Q field on the observed mi- 
crowave background in the case where the energy den- 
sity in the quintessence field is negligible during photon 
decoupling. This has been considered in ||l^]. In this 
case the effect of the dark energy is to modify the angu- 
lar distance-redshift relation . The existence of dark 
energy during the look-back time to the surface of last 
scattering shifts the acoustic peaks in the CMB power 
spectrum to smaller angular scales and larger /-values. 
The amplitude of the first acoustic peak in the power 
spectrum also increases, but not as much for quintessence 
models as for a true cosmological constant A. The basic 
features of the observed power spectrum can be fit 
pOf with either of the quintessence potentials considered 
here. Indeed, dark energy is required to reproduce both 
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the observed power spectrum and the Type la super- 
nova data. For our purposes, this look-back constraint 
is already satisfied by demanding that 17q = 0.7 at the 
present time through Eqs. (0) and (0). 



B. Energy in Q field 

There is, however, an additional possible effect of the 
Q field on the CMB which we also use as a constraint. 
If the energy density in the Q field is a significant frac- 
tion of the background energy during the epoch of pho- 
ton decoupling, it can increase the expansion rate. In- 
creasing the expansion rate can also push the / values 
for the acoustic peaks in the spectrum to larger values 
and increase their amplitude |p8| . Such an effect of an 
increased expansion rate has been considered by a num- 
ber of authors in various contexts [ p^|p^j29[| . For our 
purposes, we adopt the constraint deduced by |]l6| based 
upon the latest CMB sky maps of the Boomerang ||2^] 
and DASI collaborations that the density in the Q 
field can not exceed fig < 0.39 during the epoch of 
photon decoupling. This implies a maximum allowed 
contribution of the Q field during photon decoupling of 
Pq/pb = f^Q/(l — ^q) ^0.64. Note, however, that the 
Q field behaves differently near photon decoupling than 
during the BBN epoch (cf. Fig. [^). After the transi- 
tion from a radiation-dominated to a matter-dominated 
universe, the Q field (now coupled to the matter field) 
can contribute a much larger fraction of the background 
energy density (e.g. Curve B of Fig. |l]) than in the BBN 
epoch. Hence, the constraint on pq field from BBN must 
be considered separately from the constraint at the CMB 
epoch. 



VI. QUINTESSENTIAL INFLATION 

Another possible constraint arises if the kinetic term 
dominates at an early stage (e.g. Curve A of Fig. |l|). In 
this case pQ ~ and pq decreases with scale factor as 
. At very early times this kinetic regime can be pro- 
duced by so-called "quintessential inflation" [jllj. In this 
paradigm entropy in the matter fields comes from gravi- 
tational particle production at the end of inflation. The 
universe is presumed to exit from inflation directly into a 
kinetic-dominated quintessence regime during which the 
matter is generated. An unavoidable consequence of this 
process, however, is the generation of gravitational waves 
along with matter and the quanta of the quintessence 
field ITllJsol-lsl at the end of inflation. 



A. Energy in Quanta and Gravity Waves 

Particle creation at the end of inflation has been stud- 
ied by the standard methods of quantum field theory in 



curved space-time 
particles is just 




The energy density in created 



PB 



128 



zi 



1 



5i 



9eff{z) 



1/3 



(12) 



where Hi and zi are the expansion factor and redshift at 
the matter thermalization epoch at the end of quintessen- 
tial inflation, respectively. The factor of 128 in this ex- 
pression comes from the explicit integration of the parti- 
cle creation terms. 

When the gravitons and quanta of the Q field are 
formed at the end of inflation, one expects |jll[] the energy 
density in gravity waves to be twice the energy density 
in the Q-field quanta (because there are two graviton po- 
larization states). In this paradigm then, wherever we 
have deduced a constraint on pq, it should be taken as 
the sum of three different contributions. One is the dark 
energy from the vacuum expectation value {pq) of the Q 
field; a second is the contribution from the fluctuating 
part psq of the Q field; and a third is from the energy 
density pcw in relic gravity waves. Thus, we have 



PQ {{pq) + P&Q + PGw) 



(13) 



For the cases of interest, the energy density in gravity 
waves and quanta scales like radiation after inflation, 
Pgw-^Psq oc a~^, while the quintessence field vacuum ex- 
pectation value evolves according to Eq. (|^) which gives 
{pq) oc during the kinetic dominated epoch. This 
epoch following inflation lasts until the energy in the Q 
field falls below the energy in background radiation and 
matter, pq < ps- 

Thus, for the kinetic dominated initial conditions 
(curve A) in particular, gravity waves could be an im- 
portant contributor to the excess energy density dur- 
ing nucleosynthesis. The relative contribution of gravity 
waves and quintessence quanta compared to the back- 
ground matter fields is just given by the relative number 
of degrees of freedom. At the end of inflation, the relative 
fraction of energy density in quanta and gravity waves is 
given by 0| 



{psQ + Pgw)/pb = i/Ns 



(14) 



where Ns is the number of ordinary scalar flclds at the 
end of inflation. In the minimal supersymmctric model 
Ng = 104. Propagating this ratio to the time of nu- 
cleosynthesis requires another factor of {g-n/gi)^^^ where 
gn = 10.75 counts the number of effective relativistic 
degrees of freedom just before electron-positron annihi- 
lation, and gi counts the number of degrees of freedom 
during matter thermalization after the end of inflation. 
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In the minimal standard model gi = 106.75, but in su- 
pergravity models this increases ^ 10'^. 
Combining these factors we have 



{PSQ + PGW)/PB < 0.014 



(15) 



during nucleosynthesis. Hence, in this paradigm, the al- 
lowed values of pq/ pb consistent with nucleosynthesis 
could be reduced from a maximum of 0.056 to 0.042, fur- 
ther tightening the constraints deduced here. 



B. Gravity- Wave Spectrum 



There has been considerable recent interest 30-3^] 
in the spectrum of gravity waves produced in the 
quintessential inflation paradigm. One might expect that 
the COBE constraints on the spectrum also lead to con- 
straints on the Q field. Here we consider this possibility, 
but conclude that no significant constraint on the ini- 
tial Pq or effective potential is derived from the gravity 
wave spectrum. On the other hand, the BBN and CMB 
gravity-wave constraints discussed here can be used to 
provide useful constraints on the quintessential inflation 
epoch as we now describe. Our argument is as follows: 
The logarithmic gravity-wave energy spectrum observed 
at the present time can be defined in terms of a differen- 
tial closure fraction. 



1 dpGW 
Pc din J/ 



where the pcw is the present energy density in relic gravi- 
tons and Pc(0) = SHq/SttG — HQrn^i is the critical den- 
sity. This spectrum has been derived in several recent 
papers |^^-^. It is characterized by spikes at low and 
high frequency. The most stringent constraint at present 
derives from the COBE limit on the tensor component 
of the CMB power spectrum at low multipoles. There is 
also a weak constraint from the binary pulsar |^,^ and 
an integral constraint from nucleosynthesis as mentioned 
above. 

For our purposes, the only possible new constraint 
comes from the COBE limits on the tensor component 
of the CMB power spectrum. The soft branch in the 
gravity-wave spectrum lies in the frequency range be- 
tween the present horizon i/q = 1.1 x IQ^^^Vt^jl^h Hz and 



the decoupling frequency Vdeci^) — 1-65 x 10 "^^^m 



Hz, where we adopt JIm — 0.3 for the present matter 
closure fraction. The constraint on the spectrum can be 
written [51, 



81 

^(iH" 



gdec 
91 



1/3 / JJ^-^"^ 

m.pi 



X InM — ] < 6.9/i"2 X 10-11 



(17) 



where, fi-^ = 2.6 x 10~ /i is the present closure fraction 
in photons. The number of relativistic degrees of freedom 
at decouphng is Qdec = 3.36. As noted previously, gi is 
the number of relativistic degrees of freedom after matter 
thermalization. In the minimal standard model is gi = 
106.75. The quantity Hi is the expansion rate at the end 
of inflation. In quintessential inflation it is simply related 
to the kinetic energy pq after inflation. 



Pq(zi) = Him: 



pi 



(18) 



The logarithmic term in Eq. (17) involves present val- 
ues of the frequency Vr{0) characteristic of the start of 
radiation domination {pB — Pq a^t z — Zr), and the fre- 
quency characteristic of matter thermalization at the end 
of inflation i^i(O). The quantity vi is just, 



Similarly, 



Vr{0) = 



SO that, 



.1(0) = ^ 



Hr { Zq + 1 
Zr + 1 



Zl 



H, 



(19) 



Zl 



Zq + I 

Zr + 1 



t^r(O) 
^l(O) 



Zl + l 



(16) The identity pB — pq at z 



Ik 

128 



Hi_ 

mpi 



Zr then gives, 

2 / X 1/3 
9r 



(20) 



(21) 



(22) 



where for the cases of interest gr = 10.75. 

Collecting these terms, we can then use Eq. ( p7[ ) to 
deduce a constraint on the expansion factor Hi at the 
end of inflation 



Hi < 1.4 x lO^^TO^ 



pi 



(23) 



For kinetic dominated models, pq(zi) at the end of in- 
flation is simply related to the energy density pQ at z, 



PQiz) PQizi] 



z + 1 



(24) 



Zl + 1^ 

Similarly the background matter energy density scales as 



Psiz) = Pb{zi) 



9i 



9eff{z) 



1/3 



Z + 1 



Zl 



1 



(25) 



Considering the present energy density in photons and 
neutrinos, we can find a relation between Hi and zi: 



P7O + Pi^O 



1.1 x lO^^^^m, 



4 
pi 

1/3 



Ns if 91 1 f 
= 128^4^' 



1)" 



(26) 
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We then deduce from this and Eq. ( p3D that zi < 
8.4 X 10^^. These equations then imply that there is 
only a lower limit on pa / Pb given by the constraint on 
Hi. Combining Eqs. (pi), (pj), and (psl), we have 



Pq{z)/pb{z) 



91 



NsHl \geff{z) 



-1/3 



1 



(27) 



At our initial epoch z — 10^'^, we then deduce that , 
Pq{z)/pb{z) ^5.6 X 10^^*. This limit is not particularly 
useful because pg field must exceed at zi in order for 
the gravitational particle production paradigm to work. 
The implication is then that all initial conditions in which 
the kinetic term dominates over the background energy 
at zi are allowed in the quintessential inflation scenario. 
Hence, we conclude that the gravity-wave spectrum does 
not presently constrain the initial pq or V{Q) in the 
quintessential inflation model. 

Before leaving this discussion on quintessential infla- 
tion, however, we remark that the limits on pg / pB < 560 
derived from the BBN constraints discussed below, can 
be used to place a lower limit on the expansion rate at 
the end of inflation in this model. This in turn can be 
used to deduce a lower limit on the redshift for the end 
of quintessential inflation from Eq. (p6[). Thus, we have 



(28) 



H 

2.2 X 10^2^ < ^ < 1-45 X 10"" 
^pi 

and in the minimal supersymmetric model 



1.0 X 10-^^ < zi < 8.4 X 10' 



25 



(29) 



This nontrivial constraint then implies that the kinetic 
driven quintessential inflation must end at an energy 
scale somewhere between about 10^ and lO^'^ GeV, well 
below the Planck scale. By similar reasoning one can 
apply this argument to the gravity-wave spectrum from 
normal inflation as given in | |3^ ]. We deduce an upper 
limit to the epoch of matter thermalization oi Hi < 3.1 x 
10-^°ml, which implies zi < 7.3 x 10'^^[g{zi)/3.36]^^'^^ . 
In this case there is no lower limit from BBN as there is 
no Q field present after inflation. 



VII. RESULTS AND DISCUSSION 

The equations of motion [Eqs. (||) and (^)] were evolved 
for a variety of initial Q flcld strengths and power-law pa- 
rameters a. As initial conditions, the quintessence field 
was assumed to begin with equipartition, i.e. = 
V{Q). This implies wq = initially. This seems a 
natural and not particularly restrictive choice, since wq 
quickly evolves toward the kinetic {wq = +1) or the 
tracker solution [Eq. (pO|)] depending upon whether one 
begins above or below the tracker curve. 



Nevertheless, this initial condition does not encompass 
all possible cases. For example, one might imagine a 
somewhat strange initial condition in which Q begins 
at a very large value such that V{Q) « Pc{0) while 
Pq — is arbitrarily large. In this case, pg will 

quickly decay away as before the present epoch and 
insufficient dark energy will be present at the current 
epoch. Hence, it is possible that even though one starts 
with a large pg, one may not satisfy the present equa- 
tion of state constraint. We do not consider this as a 
likely case, however. The natural initial values for the 
quintessence field are Q < mpi, so that sufficient dark 
energy should always be present for this potential choice. 
Besides, postulating an initially small value for Q rein- 
troduces the small A problem which is what these models 
were invented to avoid. 

Constraints on a and the initial value for the Q-field 
energy density pq{z) / pb{z) at z = 10^^ were deduced nu- 
merically. These are summarized in Figure^ for both: (a) 
the bare Ratra power-law potential; and (b) its SUGRA 
corrected form. We present the quantity pq/pb as a 
more intuitive measure of the relative amount of Q-field 
energy density than to just give pq. This quantity can 
be easily converted to fig = {pq/ Pb)/{^ + Pq/pb)- For 
purposes of illustration, we have arbitrarily specified ini- 
tial conditions at z = 10^^, corresponding to T 10^^ 
K, roughly just after the time of the QCD epoch. At any 
time the energy density prei (z) in relativistic particles is 
just 



Preliz) = PjO 



3.36 

9effiz) 



1/3 



(z + iy 



(30) 



where p^o = 2.0 x 10~^^ GeV*, is the present energy 
density in microwave background photons, and we take 
geff{z) — 10.75 between z — 10^^ and the beginning 
of BBN just before electron-positron annihilation (z « 
lOiO). 

The envelope of models which obtain a tracker solu- 
tion by the epoch of nucleosynthesis are indicated by up- 
per and lower curved lines in Figures ^a and The 
general features of these constraints are as follows. If 
the initial energy density in the Q field is too large, the 
tracker solution is not reached by the time of BBN. The 
Q-field energy density can then significantly exceed the 
background energy during nucleosynthesis. This situa- 
tion corresponds to the excluded regions on the top of 
Figures ||a and 

All solutions consistent with the primordial nucleosyn- 
thesis constraints are also consistent with our adopted 
CMB Q-energy constraint as also shown at the top of 
Figures ^ and Similarly, if the initial energy in the 
Q field is too small, the universe docs not become dark- 
energy dominated by the present time. This corresponds 
to the excluded (no A) region at the bottom of the fig- 
ures. 
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The excluded regions at the top and bottom of fig- 
ures |a and |b can be easily understood analytically. 
For example, the excluded (no A) region at the bot- 
tom of these figures reflects the fact that if pq is ini- 
tially below the value presently required by f2A — 0.7 
[cf. Eq. (^)] it can not evolve toward a larger value. 
Hence, Pq/pb < 2.8 x 10"'^^ is ruled out for h = 
0.7. Similarly, the "Excluded by BBN" region comes 
from requiring that pq{zbbn) / Pb{zbbn) < 0.056. For 
this constraint we are only considering cases in which 
the Q field is approaching the tracker solution from 
above during nucleosynthesis. Hence, it is in the ki- 
netic regime in which pg cx while the background 
scales as z'^. Thus, we have pq{z = 10^^)/p_b(z = 
1012) > 0.056(10iVl0i°)'^~4 = 560 is excluded. By simi- 
lar reasoning, the "Excluded by CMB" region is given by 

Pq{zcmb)/ Pb{zcmb) > 0.64, or pq{z = 10'^'^)/pb{z = 
1012) > 0.64(10.75/3.36)1/3(1012/103)6-4 ^ 9 4 ^ ioi7. 

For the bare Ratra power-law potential (Fig. 3a) the 
main constraint is simply the requirement that the equa- 
tion of state be sufficiently negative by the present time. 
The sensitivity of the allowed power-law exponent to the 
equation of state is indicated by the wq = -0.2, -0.4, 
and -0.6 lines on Figure In the present Q-dominated 
epoch, Eq. (^0|) is no longer valid even in the tracker 
regime. Thus, the lines of constant wq must be evalu- 
ated numerically using Eq. (^. The slight slope to these 
curves comes from the fact that V{Q) / p{Q) has not yet 
reached unity, i.e. there is still some small kinetic con- 
tribution to p{Q) and the amount of kinetic contribution 
depends upon a. 

For the bare Ratra power-law potential, tracker solu- 
tions with Q!^20 are allowed if wq < —0.2. The allowed 
values for a reduce to < 9 and < 2 if the more stringent 
-0.4 and -0.6 constraints are adopted. If a is too small, 
say a < 2, then the potential parameter M becomes a 
very small fraction of the Planck mass [cf. Eq. (||)], and 
the fine tuning problem is reintroduced. 

For models in which the tracker solution is obtained by 
the time of BBN, nucleosynthesis only limits the poten- 
tial parameters if the most conservative equation of state 
limit (wq < —0.2) and most stringent nucleosynthesis 
constraint [pq/pb < 0.1%) are adapted. On the other 
hand, independently of the equation of state constraint, 
nucleosynthesis limits a large family of possible kinetic- 
dominated solutions in which the Q-field energy density 
exceeds the background energy prior to and during the 
nucleosynthesis epoch. Such models are excluded even 
though they provide sufficient dark energy in the present 
epoch. 

For the SUGRA-corrected Q fields (Fig. 3b), the con- 
straint from Wq is greatly relaxed. In fact, wq is suffi- 
ciently negative {wq < —0.6) for all a < 10''. The reason 
is that all tracker solutions have wq « — 1. This is be- 
cause Wq decays much faster toward -1 for the SUGRA 
potential. Also, the potential has a finite minimum which 



is equal to the present dark-energy density. The Q field 
quickly evolves to near the potential minimum and has 
negligible kinetic energy by the present time. Any po- 
tential which becomes flat at late times gives wq sa — 1 
and the dark energy looks like a cosmological constant 
A. All SUGRA models which achieve the tracker solu- 
tion also have a small pg during primordial nucleosyn- 
thesis. Hence, there is no constraint from nucleosynthesis 
except for those kinetic-dominated models in which the 
Q field is still far above the tracker solution during the 
nucleosynthesis epoch. 

We do note, however, that if a lower limit of wq > — 0.8 
is adopted for tracker solutions from then only a 
power law with a^30 is allowed as indicated on Figure 
3b. This makes the SUGRA potential the preferred can- 
didate for quintessence. The large a implies values of 
M in Eq. (||) close to the Planck mass, thus avoiding 
any fine tuning problem. However, this quintessence po- 
tential is constrainable by BBN should the light-element 
constraints become sufficiently precise to limit pg at the 
0.1% level. 

VIII. CONCLUSIONS 

We conclude that for both the bare Rata inverse power- 
law potential and its SUGRA-corrected form, the main 
constraints for models which achieve the tracker solu- 
tion by the nucleosynthesis epoch is from the require- 
ment that the equation of state parameter becomes suf- 
ficiently negative by the present epoch. The main con- 
straint from nucleosynthesis is for models in which the 
Q field has not yet obtained the tracker solution by the 
time of nucleosynthesis. The SUGRA-corrected potential 
is the least constrained and avoids the fine tuning prob- 
lem for M . Therefore, it may be the preferred candidate 
for the quintessence field, although BBN may eventually 
limit this possibility. We also note that the constraints 
considered here provide useful constraints on the regime 
of matter creation at the end of quintessential inflation. 
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FIG. 1. Examples of the evolution of energy density in pQ 
and the background fields ps as a function of redshift for 
an inverse power-law effective potential with a = 5. The 
BBN epoch is indicated by a short line segment. The location 
of the transition from radiation to matter domination (for 
Q,M ~ 0.3), and the CMB epoch, are also indicated. Curves 
A and C illustrate cases in which the Q field does not achieve 
the tracker solution until well after the BBN epoch. Curve A 
is excluded, curve C is not. Curve B is an example of a tracker 
solution which is allowed by the nucleosynthesis constraints. 
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FIG. 2. Constraints on the ratio of the energy density 
in the quintessence field to the background energy density 
pq/pB (at T = 1 MoV) from the primordial abundances as 
indicated. The allowed region corresponds to 1^ < 0.247 and 
D/H < 4.0 X 10"^ 



10 



(a) 



(b) 




a 

FIG. 3. Contours of allowed values for a and initial pq/ pB (at z — 10^^) from various constraints as indicated for (a) the 
bare power-law potential, and (b) the SUGRA corrected potential. Models in which the tracker solution is obtained by the 
BBN epoch are indicated by the upper and lower curves. Values of a to the right of the lines labeled wq — —0.6, —0.4, —0.2 
on (a) are excluded by the requirement that the present equation of state be sufficiently negative. The BBN constraint for a 
maximum energy density in the Q field of 0.1% (dotted line) and 5.6% (solid line) are also indicated. For the SUGRA potential 
(b) all tracker solutions to the right of the region labeled wq = —.8 are allowed. 
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